Godunov mixed methods on triangular grids for advection–dispersion equations

A time-splitting approach for advection–dispersion equations is considered. The dispersive and advective fluxes are split into two separate partial differential equations (PDEs), one containing the dispersive term and the other one the advective term. On triangular elements a triangle-based high resolution Finite Volume (FV) scheme for advection is combined with a Mixed Hybrid Finite Element (MHFE) technique to solve dispersion. This approach introduces an error proportional to the time step and the overall scheme is only first order accurate if special care is not taken in the definition of the numerical flux approximation for advection. By incorporating the diffusive effects into the definition of this numerical flux, near second order accuracy (up to a logh factor) can be proved theoretically and validated by numerical experiments in both oneand two-dimensional cases.